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ABSTRACT 

The Laser Interferometer Space Antenna {LISA) will detect gravitational wave 
signals from coalescing pairs of massive black holes in the total mass range (10^- 
10^)/(1 +z) Mq out to cosmological distances. Identifying and monitoring the elec- 
tromagnetic counterparts of these events would enable cosmological studies and 
offer new probes of gas phys ics around well-characterized massive black holes. 



Milosavljevic & Phinneyl (|2005|) proposed that a circumbinary disk around a binary 
of mass ~ 10^ Mq will emit an accretion-powered X-ray afterglow approximately one 
decade after the gravitational wave event. We revisit this scenario by using Green's 
function solutions to calculate the temporal viscous evolution and the corresponding 
electromagnetic signature of the circumbinary disk. Our calculations suggest that an 
electromagnetic counterpart may become observable as a rapidly brightening source 
soon after the merger, i.e. several years earlier than previously thought. The after- 
glow can reach super-Eddington luminosities without violating the local Eddington 
flux limit. It is emitted in the soft X-ray by the innermost circumbinary disk, but it 
may be partially reprocessed at optical and infrared frequencies. We also find that the 
spreading disk becomes increasingly geometrically thick close to the central object 
as it evolves, indicating that the innermost flow could become advective and radia- 
tively inefficient, and generate a powerful outflow. We conclude that the mergers of 
massive black holes detected by LISA offer unique opportunities for monitoring on 
humanly tractable timescales the viscous evolution of accretion flows and the emer- 
gence of outflows around massive black holes with precisely known masses, spins and 
orientations. 

Subject headings: accretion, accretion disks — black hole physics — gravitational 
waves — quasars: general 



1. Introduction 



Massive black holes (MBHs) are abundant at all observabl e redshifts, manifesting them- 
selves most often as active galactic nuclei (AGN) and quasars (|Kormendy & Richstond Il995 : 
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Richstone et al.lll998|) . A me rger of a pair of galax ies each containing a MBH will result in the 
formation of a MBH binary (IBegelman et al.lll980|) . Given the hierarchical structure formation 
histo ry of the universe, and that most or all formed galaxies appear to harbor a MBH in their nuclei 
(e.g., ^Magor rian et al. 1998k there should be many MBH binaries formed throughout cosmic time 
(lHaehnelt.il994l:lM enou et al.l200ll : lwonteri et al Jbooilwvithe & LoeblboOsUsesana et al.ll2007l : 



Tanaka & Haiman.200 9). Once compact, these binaries rapidly lose orbital energy via gravitational 
radiation (GW) and coalesce. The Laser Interferometer Space Antenna {LISA) will detect gravita- 
tional waves resulting f rom MBH me rgers of binaries with total mass M = ( 10^ - 10^)/( 1 +z)Mq out 
to redshift z ~ 20 (e.g.. ISchutal2009() . LISA is expected to be able to constrain the individual red- 
shifted mass {l+z)M, the mass ratio ^ < 1, spins of the MBHs, and the luminosity distance of the 
source with precision, thereby providing an unprece dented test of general relativity and probing the 
assembly history and demography of MBHs (e.g., lHughesll2002l : IVecchioll2004l : Lang & Hughes 
20061) . However, the space-based detector will be unable to measure the redshift (due to the fun- 
damental degeneracy between source mass and redshift) or the precise angular location (due to the 
fact that sources are located using the mod ulation of the signal due to orbital motion and se para- 
tion of the detector elements ) of the source Jcutleill 19981 : IhoIz & Hugheslboosi: IkocsIs et al.ll2006. 
2007l : lLang & Hughesll2008l) . 



The coalescence of a MBH binary is not by itself expected to emit an observable electromag- 
netic (EM) signal. If, however, MBH mergers have associated and readily identifiable EM signa- 
tures, their detection would comp l ement LISA by helping determine soiirce redshifts an d locations 
on the sky (IhoIz & Hugheslboosi : IkocsIs et allliooel : bloom et~aDb009l : lphinnevlb009l) . An obvi- 
ous candidate to effect such emission is the abundant gas linked to galactic MBHs. Rapid, sustained 
gas accretion onto MBHs is required to explain the luminosities and number of observed AGN and 
quasars. Numerous studies suggest th at a gas -rich environment aids the formation of close MBH 
binaries in recently merged galaxies (lEscala et al.l |2005|; lOotti et al.l 120071 : ICallegari et al.l l2009l : 
Colpi et al.ll2009i) . and may facilitate the further shrinking of the orbit to where GW emission can 
enact the merger (Ivanoy et al. [l999l: Armitage & Nataraian 2002; MacFadyen & Milosavljevic 



20081 : ICuadra etal.lb009l : 



Lodato et al 



20091) . Various types of interactions between MBH bi- 



naries and their gas environments have been proposed as viable mechanisms for observable EM 
emission. Numerous studies have calculat ed the EM e r nission frorn the r esponse of the gas disk 
to the mass loss and gravitatio nal recoil (PeresI 19621: Bekenstein 1973) effects which accom- 
pany the merger process (e.g.. Bode & Phin nev boC Tzl; Schnittman & KrolikI 2008: Lippai et al. 



2008 



Shields & Bonning 



Rossi et al.„2009l : 



2008 



O'Neill et all 120091: 



Anderson et al. 



20091). IChang et al. 



Megevand et al.li2009i : iCorrales et al.ll2009l : 



(|2009h have suggested a lum inous EM sig 



nal from tidal and viscous excitation of fossil gas trapped inside the binary's orbit. I KrolikI (120101) 
pointed out that even small amounts of gas that is present in the immedi ate vicinity of a merg- 
ing binary can power short afterglows on Eddington luminosity scales. iKocsis & Loeb (I2008.) 
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have proposed an infrared a f terglow from the dissipa tio n of GW through a surrounding gas disk 
(see, however, lKrolikll2010|) . IPalenzuela et alj (|2009alJbh : lM6sta et alj (|2009|) have emphasized the 
possibility of variable emission due to the perturbation and enhancement of the local EM fields. 
Haiman et alj (l2009a f) have raised the possibility of detecting the binary as a periodic variable 
source prior to the merger. 



In this paper we revisit the a ccretion afterglow mechan ism proposed by lMilosavljevic & Phinney 
(2005; hereafter MP05; see also Liu et al. 2003; Liu 2004) for a geometrically thin circumbinary 
disk. We summarize their model as follows: 



The tidal torques from the binary open a gap in the circumbinary gas. The gas inside the 
gap accretes, while the gas outside is held at bay by tidal torques. What results is a cir- 
cumbinary disk with the binary inside the central cavity. Because it is largely "missing" the 
inner, highest-temperature region, the disk ca nnot easily produce thermal X-rays via viscous 
dissipation (see, however, Chang et ahl 2009). The system can remain in this configu ration 
for an extended period of time ( Ivanov et al.lll999 ; MacFadyen & Milosavljevi^boosb as the 
timescale for the extraction of angular momentum from the binary's orbit is long. 

Once the binary reaches a semimajor axis a ~ lOOGM/c^, it rapidly loses orbital angular mo- 
mentum and energy via G W emission and the orbit closes faster than the surrou n ding gas can 



visco usly follow (see also lArmitage & NatarajanI |2005| ; iHaiman et al.l l2009al ; IChang et al 



20091). The binary coalesces, producing the GW signature detectable by LISA. 



3. The gas (no longer held back by binary torques) reaches the center of the cavity by viscously 
spreading, and the corresponding accretion flow, deep in the potential well of the MBH 
remnant, emits an X-ray afterglow. 



By taking the difference between the time it takes the binary to merge and the time it takes for 
the bulk of the gas at the cavity edge to reach the central remnant, ImPOSI estimated that the X-ray 
afterglow would occur ?em ~ 7(M/1O^M0)^-^(1 +z) yr after the GW signal (in the observer's rest 
frame). 

In this study, we use the term "accretion afterglow" to denote the emission described above. 
We stress that this is a distinct mechanism from the various afterglow mechanisms that are powered 
by the response of the accretion disk to the mass loss and/or recoil of the central MBH remnant 
(references above). In the latter scenarios, mass loss and/or recoil introduce additional eccen- 
tricities in the gas orbits, and the emission is powered by the circulariza tion of the orbits and 
shock/wave dissipation. In the accretion afterglow scenario considered by ImPOSI and in this pa- 
per, the emission results from the deepest parts of the MBH potential well becoming accessible to 
the circumbinary disk as the binary shrinks and merges. Before the merger the disk is deficient 
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in the hard UV and X-ray frequency range because of the central cavity, while after the merger 
this high-frequency emission emerges as the cavity fills with gas. Although distinct, the various 
classes of emission mechanisms considered are not entirely unrelated. For example, mass loss 
introduces a roughly constant eccentricity everywhere in the disk, and the resulting luminosity of 
a circularization-powered afterglow (references above) would dep end delicately on the innermost 



density profile of the circumbinary disk at the time of merger (see lCorrales et al.ll2009 l for a com- 
parison of afterglows for different density profiles). A high surface density for the innermost disk 
(where the avai lable specific orbital energy is highest) at merger would enhance the mechanism 



emphasized by lO'Neill et al.l (120091) while it would suppress the pre- and post-merger emission 
contrast considered in this paper; a low surface density would have the reverse consequences for 
the two mechanisms. 

The purpose of this work is to investigate the onset and time dependence of the afterglow 
in greater detail. We use an idealized semi-analytic framework to model the time dependence of 
the electromagnetic signal. We present explicit, integral solutions for the time evolution of one- 
dimensional (viz . geometrically thin and azimuthally symmetric), viscous Keplerian disks. (See 



similar work by iShapird 120101 ) Our calculation method is valid insofar as the disk kinematic 
viscosity can be described as a simple function of radius, and as long as after decoupling the 
inner gas is minimally affected by the gravitational torques from the binary. This last point is 
worth emphasizing, as in general the decoupling of the binary from the gaseous influence does 
not guarantee that the gas is entirely free from the binary influence. Even after GW emission 
has become the dominant mechanism driving the orbital evolution of the binary, the binary's tidal 
torques can still influence any gas that is able to remain in the vicinity of the binary's orbit. 

This paper is organized as follows. In §(2] we review the disk properties as the GW-driven 
closing of the binary decouples from the viscosity-driven spreading of the circumbinary disk, and 
introduce a semi-analytic model for the subsequent evolution of the disk's surface density. Various 
derivations and intermediate results used in this section are provided in the Appendices |A]& [HI We 
present in §|3]the predictions of the model for the light curve and spectral evolution of the resulting 
electromagnetic emission. We also discuss in that section the possibility that the X-ray afterglow 
could be reprocessed by the outer disk, as well as possible effects of advection as the inner disk 
becomes geometrically thick. We conclude in ^ 



2. Modeling the Binary-Disk System 



Below, we first discuss the properties introduced by iMPOSi for the evacuated circumbinary 
gas disk at the time when the evolution of the binary orbit becomes dominated by GW emission 
and thus decoupled from the evolution of the surrounding disk, which evolves viscously (see also 
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Haiman et alj|2009b|) . We then describe the semianalytic integral formalism for the viscous evolu- 



tion of an accretion disk with an arbitrary initial surface density profile. Additional derivations and 
intermediate results for both the disk properties and the semianalytic model are detailed in Appen- 
dices A and B, respectively. Throughout this paper, R is the radial distance from the center of mass 
of the binary in the plane of its orbit; a is the binary's semimajor axis; M is the total mass of the 
binary; and q < I is the binary mass ratio. The quantities c, G, nip, h, k, a and ctt are the speed of 
light, the gravitational constant, the mass of the proton, Planck's constant, the Boltzmann constant, 
the Stefan-Boltzmann constant and the Thomson cross section, respectively. Times, frequencies 
and rates are in the rest frame of the binary, unless noted otherwise. 



2.1. The Circumbinary Disk at Decoupling 



The "cold" gas in the nucleus is assumed to settle into a geometrically thin, rotationally sup- 
ported disk. In this paper, we are concerned with the properties of the innermost disk and the binary 
in the very last stages of its evolution. We assume that the angular momentum of the d isk is aligned 
with the orbital angular momentum of the binary (i.e., no warps; llvanov et al.lll999l ; however we 
shall revisit the possibility of warping later), and that the binary has been circularized by inter- 
action s with the surrounding gas (see lArmitage & NatarajanI llOOSL iMacFadyen & Milosavljevic 
20081 and lCuadra et al.ll2009l for possible complications). 



(2.1) 



The disk is able to evolve locally on length scales ~ on the viscous timescale 

3 Z/ 

where iy{R) = ^aPgas/ip^) = ^akT / (ixnipQ) and PgasjP, T and fi = 0.6 are the pressure, density, 
temperature and mean molecular weight of the gas, respectively. The quantity is the angular ve- 
locity, assumed to be approxim ately Keplerian: Q ^ y/G M/R^. The parame ter a is the viscosity 
parameter in the disk model of ( Shakura & Sunyaev 19731) . We foUow lMPOS in adopting a viscos- 
ity prescription where u scales only with the gas pressure, as opposed to the total (gas, radiation 
plus magnetic) pressure. This choice is motivated by calculatio ns that suggest that disks where vis 



cosity scales with t he total pressure may be thermally unstable (jShakura & Sunyaevll 19761: IPringle 



19761 : see, however. iHirose et al.ll2009 



lay De 

3.0 



At a distance of R^ 2a, a gap is opened in the circumbinary d isk through tidal interact ions 
(see, e.g. JArtymowicz et al. I ll99l[ IMacFadyen & Milosavljevii2008l and lHaiman et al ] l2009bh . In 



' However, radiation p ressure dominated disks may still be viscously unstable (ILightman & Eardley 



1974 



Pira 



2 



1978[lHirose et al.ll2009al) . Because the afterglow mechanism discussed in this paper takes place in the radiation- 
dominated part of the disk, the findings reported here are contingent on the answer to this open theoretical question. 
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this region the input of angular momentum from the tidal torques exerted by the binary counteracts 
the viscous spreading of the disk. F or simplicity^ we n eglect the small amount of fossil gas that may 
be present inside the binary's orbit (|Chang et al.l2009() . The gap, at first annular, becomes a roughly 
circular central cavity after the inner gas is consumed. The boundary between the cavity and the 
disk is characterized by a steep surface density gradient, a pile-up of gas caused by tidal decretion 
in the vicinity of the binary's orbit and viscous accretion outside the decreting region. Whereas a 
steady-state accretion disk has a mass accretion rate M{R) = Sttz/S that is constant with radius, a 
decretion disk with a steady inner torque instead satisfi es constant viscous torque, i.e. 37ri?^z/SQ = 
constant (e.g., iRudak & Paczvnskilll98ll : IPringlelll99ll ). It is a priori unclear how close to accretion 
or decretion a certain syster n is at decoupling, as this will depend on the system parameters and its 
past accretion history (e.g. Jlvanov et al.ll 19991 : IChang et al.ll2009() . We thus consider various pos- 
sibilities for the circumbinary disk density profile at decoupling. Some of the innermost gas may 
be still able to accrete across the binary's orbit in narrow streams, but the accretion rate is expected 
to be only a small fraction of that of a comparable accretion disk around a single black hole (e.g.. 



20081) 



Lubowetal.lll999l : iGunther et al.ll2004 iHavasaki etal.ll2007L l2008l : iMacFadven & Milosavljevic 



The timescale for the binary separation to shrink via GW emission is 



da/dt 16 C 



4.90 yrx 



lOOGM/c 



(2.2) 



where C,=Aq/{l+qfh the symmetric mass ratio scaled to unity for equal-mass binaries. Equation 
(12.21) is valid for binaries with non-extreme mass ratios and circular orbits. Because fcw ^ ^visc 
initially in the circumbinary disk, the gas is able to respond promptly to the gradual and relatively 
slow shrinking of the binary, probably maintaining a roughly constant geometric ratio between the 
disk inner edge and the orbital semimajor axis, i.e. A =R/{2a) ~ 1. 

Because the disk viscosity z/ is a weak function of radius outside the inner edge (see Appendix 
lAl) . roughly speaking the viscous time there scales as ?visc(^) oc R^. As the binary closes, tQ^{a) oc 
will inevitably become shorter than the disk viscous time at the inner edge, and the binary will 
begin to close faster than the bulk of the gas at the inner edge can follow. The subsequent evolution 
of the binary is driven by GW emission and thus causally decoupled from that of the disk. In the 



absence of binary torques, the gas at the cavity edge wou 
the approximate boost factor /3 ~ 0. 1 wa s introduced by 



d fill th e cavity in a time ~ /3?visc, where 
t o acco unt for the limit of a very 



MP05 



steep surface density gradient at the edge ( Lynden-Bell & Pringle 1974 1). Decoupling occurs, then, 
when tcwif^i) ~ /3?visc(2Aa). This condition, along with the assumption that the viscously dissipated 
energy is locally emit ted as t hermal radiation, specifies the various properties of the inner edge of 
the disk described by ImPOSI . We closely reproduce their results in Table [U and detail the various 
intermediate calculations in Appendix |Al Note that while we recover the parameter dependencies 
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exactly, the magnitudes of our disk properties differ somewhat from those of lMPOSl This is due to 



differences in how we hav e calcu lated the thermal structure of the disk, the most significant being 



that the flux expression in IMPOSI is a factor of two lower than ours (compare their equation [2] 
and preceding formula for with our equation [ (13.41) 1 below). Throughout the rest of this paper, 
we use the subscript "0" to denote the value of a quantity at decoupling. When applied to a disk 
quantity with radial dependence, the subscript shall also denote the value at the inner edge of the 
disk. For example, is the semimajor axis of the binary orbit at decoupling and is the angular 
velocity at the inner edge of the disk at decoupling. 

The surface density at the inner edge of the circumbinary disk at decoupling cannot be de- 
termined through scaling arguments alone. However, we expect it to be greater than the value 
expected for a steady thin accretion disk. Just outside the edge, the gas may resemble a decretion 
disk, with Bttz/S oc R~^I^ approximately, because of mass accumulation. Sufficiently far from the 
binary and its tidal torques, we expect the gas to behave more like a steady thin accretion disk, with 
a local mass supply rate M = dTivT, r oughly constant and determined by external conditions (e.g., 
Ivanov et aDll999 : Chang et al. I I2OO9I) . It follows that if StivT, is monotonic outside the edge then it 



must be greater than the outer mass supply rate of the disk. Therefore, we parametrize the surface 
density at the inner edge via the arbitrary relation So = 5A/Edd/(37rz/o), where 5 is a dimensionless 
parameter and MEdd = A-nGMm^/ (rjcaj-) is the Eddington accretion rate onto a MBH of mass M for 
a radiative efficiency 77 = 0. 1 . Effectively, S is the product of the gas supply rate onto the disk in 
Eddington units and the enhancement of the local surface density at the disk edge due to the binary 
torques and mass accumulation^ In general, S can exceed unity even if the disk mass supply rate 
and the local accretion rate are both sub-Eddington. 



2.2. A Simple Model for the Viscous Evolution 

Followin g deco upling, the binary outruns the surrounding gas to the center and merges in a 



time ?Gw,o/4. iMPOSi estimated that an accretion-powered X-ray afterglow would take place when 
the central cavity is filled, viz. after a time t^M = /3^visc,o-^Gw,o/4 ~ (3/4)?gw,o after the GW event. 
Below, we revisit this estimate by calculating the time dependence of the signal with a simple 
model for the disk evolution. We proceed by assuming that following decoupling the binary's orbit 
closes faster than the gas can follow, viz. viscosity is the only source of torque in the disk and 



^ For practical purposes S is equivalent to the parameter "m" used in lMPOSl We avoid using the latter notation to 
prevent confusion with the actual local accretion rate, which is discussed below in ^2.21 Also note that So and S do not 
relate linearly (see Table [T]l- The relationship is more complex because the disk gas viscosity i/q and the decoupling 
orbital radius oq also depend on S through the midplane temperature. 
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the g a s orbits remain circular and Keplerian. Numerical simulations (e.g.. lArmitage & Natarajan 
2OO5I : iHavasaki et all I2OO7I : iMacFadven & Milosavljevi^boOgl : Icuadra et allboool) of MBH cir- 
cumbinary disks indicate that both the eccentricity and the deviation of the angular velocity from 
the Keplerian value are small (of order a few percent) at the radius where the binary tides truncate 
the disk. Likewise, the mass reduction of the binary at merger due to GW emi ssion will only intro 



duce orbital eccentricities (roughly equal to the f ractional mass loss) of > 10 ^ (|Tichy & Marronetti 



20081 : lO'Neill et al.ll2009l : iReisswig et al.l 120091) . At radii and binary masses of interest here, or- 
bital speeds are too great for the gas to be significantly affected by gravitational recoil of the 
MBH remnant: Q ~ ^jRc^/GMc ^ Vrecoii ~ 300 km s~\ After the binary has merged, the central 
potential is that of a rotating black hole and circular orbits near the center will have somewhat 
super- Keple rian angular velocities . Approximate "pseudo-Newtonian" expressions of the Kerr po- 
tential (e.g., Artemova et al. 1994 ) suggest that for a merger remnant with less than maximal spin, 
the deviation from the Keplerian value is of order ten percent at the radius of innermost stable 
circular orbit and falls off as roughly oc . All of these modest deviations and perturbations are 
neglected in our models. 

The viscous evolution of a geometricall y thin, cylind rica lly symmetric Kep lerian disk can be 
described by the standard equation (e.g. see iPringlelll98ll and Frank et al.ll2002h 



d d 

2'ITR-m:t) = — 

ot oR 



dR ^ ' 



(2.3) 



where S is the surface density of the disk and the innermost derivative on the right-hand side 
describes the angular momentum gradient. The left-hand side of (12.31) describes the rate of change 
of the radial mass distribution, while the right-hand side gives the radial gradient of the local mass 
flow 

M = 37rz/S(l+2m), (2.4) 

where m(r) = d\n{ijTj)/ dXwR. In the above equation, the dimensionless factor in parentheses 
equals unity for steady-state accretion solutions (m ^ 0) while it vanishes in steadily decreting 
solutions (m ^ -1/2). 



If the kinematic viscosity is a function of radius only, equation (12.31) is a linear diffusion 
equation. In the special cas e u oc i?", the solution for subsequent viscous evolution (see, e.g. 
Lynden-Bell & Pringldl 19741 : a derivation is provided in our Appendix |B]) is given by 



poo 

T.{r,t)= j G{r,r\t)i:{r\t = Q)dr', 
Jo 



(2.5) 
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where 



G(r,/,0 



2 — n 



^-l/4-«/5/4/. 



X exp - 



l/(4-2n) 



2/l-«/2^1-«/2 



l-^/-exp 



2-n 



(2.6) 



T^i.T^ ''"merge ) _ 

is the Green function, S(r, ? = 0) is the initial surface density profile (at decoupling), is the 
modified Bessel function of the first kind, TZ is the ramp function, and r = R/Rq, r' = R' /Rq and 
r = 8(1 -n/2)^?/?visc.o are dimensionless variables. At decoupling, r = 0; the binary merges when 
2/3(1 -n/2)^; roughly speaking, gas co ncentrated at r = 1 diffuses to the center in a time 



T = Tr 



merge ■ 



r ~ 1 (hence the value /? ~ 0. 1 adopted by lMP05|) 



The second term in the curled brackets of equation (12.61) accounts for the fact that, after 
merger, the MBH remnant imposes a zero-torque boundary condition at some finite radius R^ = 
Ror^,, which we associate with the remnant's innermost radius for marginally stable circular orbits[^ 
(RKca). This ra dius depends on the MBH spin as a.^i,, = (1 /3)^yRiscoC^ /GM(4- y/3RiscoC^/GM-2) 
(|Bardeenlll970l) . where -1 < aspin < 1 is the standard dimensionless spin parameter of the MBH. 



For test particle orbits, one finds R^ = 6GM/c^ if agpm = 0, and 7?* = GM/c^ if agpm = 1 and the spin 
of the hole is perfectly aligned with the test particle's orbital angular momentum. Immediately 
after the merger, the MBH remnant is expected to h ave moderately high but non-maximal spin for 



a wide range of physica lly plausible scenarios (e.g. iHughes & BlandfordI 12003 



Berti & Volonteri 2008 



Baker et al. 



King & Pringle 



200 



i 



20061 



see however the "chaotic" accretion scenario of 
which predicts lower MBH spins). We thus adopt the intermediate value i?* = 3GM/c^, which is 
accurate to within ~ 20% in the range of spin 0.65 ^ aspin ^ 0.9. Our solutions given by equation 
(12.61) depend only on the initial surface density profile (which may be arbitrarily complex and need 
not be differentiable), as well as on the values of n < 2 and r*. 



Figure [T] shows time-dependent solutions calculated following equation (12.51) for two differ- 
ent initial surface density profiles at decoupling (upper left panel). In solid lines, we evolve for 
demonstrative purposes the plainest possible model: a surface density profile obeying the decre- 
tion power-law S oc r"^/^"", truncated by a step function at the cavity edge. We also evolve a 
second, more plausible profile (dotted lines): 

E(r,0) 1 



r«+V2 + exp[-11.9A(r-l)]' 



(2.7) 



This a nalytic profile is inspired by the hydrodynamic simulation results of lMacFadyen & Milosavljevic 
(l2008h . who found an azimuthally averaged inner surface density profile with S(i? ^ 2a) oc exp(-5.95i?/a) 



^See, however. iKrolik & Hawlevl(l2002b . who showed that there are several ways to define an effective inner radius, 
some of which are quite different from /?isco- 
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for a thin disk around an equal-mass binary. Equation (12.71) reproduces this exponential behavior 
for r ^ 1 and the decretion power-law behavior S ^ Sor"^/^"" for r > 1 . Note that this profile does 
not peak at the inner edge r = 1, nor is its peak value equal to the parametrized quantity Sq. In 
the context of equation (12. 7L the value of Sq sets the scale for the power-law region of the surface 
density profile just outside the inner edge. 

In Figured! the physical length scale and the time with respect to the binary merger is calcu- 
lated using the relations in Table[I]and the fiducial parameter values = q = a_i = A = 5 = 6*0.2 = 
1. We take as n = 0.4 as our fiducial viscosity power-law. For our disk solutions, this value leads 



to agreement between the late-time disk profile vTi oc \- ^JR^./R and the physical viscosity pre- 
scription u oc Pgas/ip^) of the ftgas disk. (We refer the reader to Appendix lAl for a brief discussion 
on the parameter n.) In Figure [2l we also evolve a second case with n = 11/170 ~ 0.065, which 
is consistent with v oc PgaJ(p^) at early times just outside the inner edge for a disk that has the 
decretion profile -1/2 there. In both figures, we graph the profiles at decoupling (3.1 years 
before the merger), 1 year before merger, at merger, and 9.2 years af ter me rger. The last snapshot 



corresponds to a time /3fvisc,o after decoupling, our re-evaluation of lMPOSF s nominal estimate for 
the onset of the X-ray afterglow. 

These two figures suggest that the qualitative evolution of the disk does not depend very 
sensitively on n. The viscous evolution timescale governing the Green's function scales as (1 - 
and thus the temporal dependence on n is small as long as the quantity 1 -n/2 is of order 
unity. In both figures, the qualitative evolution is such that at early times and large radii {t < 
tyisc(R)), the disk maintains the initial mass distribution. At late times and small radii (t > fvisc(^)), 
the disk approaches an analytic quasi-steady density profile. Before the merger, we do not specify 
an inner boundary condition and so the quasi-steady surface density profile satisfies E oc z/"' oc R~". 
After the merger, however, we assume that the MBH remnant imposes a zero-torque boundary 
condition at finite radius R^, which leads to a quasi-steady surface density profile satisfying S oc 
(1 - ^/R^JR)R~". The two initial profiles evolve to become largely indistinguishable even before 
the binary merges. For a variety of initial profiles, we have confirmed that the solutions converge 
on timescales shorter than the time to merger if the profiles share the same value of A, which is 
closely related to the location of the maximum of the surface density gradient, and have identical 
profiles just outside the inner edge. 

We note that the physical circumbinary disk is not likely maintain a m -1/2 decretion profile 
at radii far outside of the inner edge. Where the gas is causally detached from the binary torques 
and the mass accumulation of the inner edge, it is expected to have th e m ^ profile of a steady 



accretion disk. The simulations of iMacFadyen & Milosavljevid (|2008r) suggest that this transition 



to a steady accreting solution does not occur until the radius is at least several times the size of the 
central cavity. Because the viscous time scales as ?visc(^) oc R^~'\ we do not expect the gas profile 
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at decoupling at these outer radii to contribute to the EM afterglow on timescales of observational 
interest. 

Despite decoupling, some gas is able to follow the binary closely. This raises the possibility 
that some gas could continue to partially influence the subsequent orbital evolution of the binary. 
Such a dynamical effect will be proportional to the gaseous mass present in the vicinity of the 
binary. Since the surface density profile at decoupling cannot viscously evolve significantly at radii 
beyond several times Rq, we choose to conservatively evaluate the maximum dynamical influence 
of the circumbinary gas based on the mass inside IORq. All of the disk profiles in Figures [T] and [2] 
have a total mass inside R < IORq of 

Mdisk(i? < 10i?o) ~ 10"^M X a-_\-^^S'^-°^X^-^Ml-^\PCf'^0-^-^\ (2.8) 

where the parameter dependencies above are calculated from Table \T\ and Appendix |A] using the 
fact that Mdisk oc EqA^Qq. The low, conservative value of the disk to binary mass ratio in equa- 
tion (12.81 ) suggests that the disk is unlikely to have a dynamical effect on the binary merger, 
and to contaminate the GW signal at levels detectable by LISA. (Note, however, the scenario of 



Armitage & Natarajanll2005L in which the disk can imprint L/5A-observable perturbations on the 



binary's orbit prior to the decoupling time.) However, such gas may stil l contribute an observable 



EM s ignature because the available specific gravitational energy is large (|Chang et al.ll2009l : lKrolik 



2010h 



Even if the gas cannot influence the binary, it is possible for the binary to continue to influence 
the evolution of the innermost gaseous disk, even after decoupling. Our models assume that the 
evolution of the gas and that of the binary can effectively be treated independently, i.e. that the 
region where the potential is non-Keplerian shrinks faster than the gas can follow. Although this 
is a useful first approximation, it is not strictly guaranteed by a decoupling criterion based on GW 
and viscous timescales, as defined earlier. Indeed, our solutions allow for gas as dense as > 10"^ So 
to flow to R^ a{t) prior to the merger (see, e.g., Figure[IJ)). To understand why this happens, one 
can compare the radial gas velocity Vr = -M / (InRT,) = 3{l + 2m)i'/{2R) with the shrinking rate 
of the hypothetical disk edge, d{2Xa)/dt. The decoupling condition uq / {da / dt)^ = /3fvisc,o directly 
implies d{2Xa)/dt = /3~'7?o/?visc.o = 3/3"Vo/(2i?()). We thus find that the nominal radial velocity 
(ignoring binary torques) at the inner edge at decoupling is given by the following expression: 

d(2Xa) 

y«.o = 0.1/3_i(l + 2mo)^— ^. (2.9) 

dt 

Above, niQ is the index (91n(z/S)/(91ni? evaluated at the disk edge, which in principle may be 
much greater than unity. For sufficiently steep surface density profiles, the innermost gas can 
continue to flow into the binary's vicinitjoeven after decoupling as defined earlier. Thus, the binary 



*In fact, because our viscosity prescription is not causal, Vr can be supersonic for large positive values of niQ. Flow 
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can continue to open a gap in the disk for some time until the binary is truly able to outrun the 
surrounding gas. This happens when the binary's orbital radius has shrunk below the decoupling 
value ao by a factor [0. 1 /3_i(l + 2m)]^/^"'^^\ The profiles given by the analytic form in equation (12.71 ) 
have mo > 10 for each of the viscosity prescriptions shown in Figures [T] and |2l We thus estimate 
the binary can outrun the gas in these profiles when a ^ O.Saoj3j/^"^^\ at a time approximately 
IMg^ yr before the merger. This suggests that the time when the binary torques cease to influence 
the circumbinary gas may be as late as ~ 0.05?visc,o before the merger. 

To evaluate the possible effect of the binary's continued influence on the gas after formal 
decoupling, we consider a scenario where the inner gas has maintained an exponential inner density 
profile similar to the one in equation (12.71) . scaled to the new semimajor axis of the binary, until one 
year before the merger as estimated above. This revised configuration is shown as dashed lines in 
Figures [Hand [21 It is evolved viscously in the same way as the other two scenarios in panels (c) and 
(d) of each figure. While the first two classes of profiles represent the assumption that the binary 
ceased to tidally interact with the circumbinary gas at decoupling (tew = PUisc), the third class of 
profiles is evolved assuming the tidal interaction continues for a while longer after decoupling, 
until the binary is able to truly outrun the disk edge as estimated above. We find that the three 
types of profiles are unlikely to be observationally distinguishable from each other (compare long- 
dashed and dotted lines in Panels b, c, d in Figures [T]and[2l); there are small differences until merger 
time, but they are quickly wiped clean by the boundary condition associated with the single MBH 
remnant. Figures [I]and[2]illustrate that neither the density gradient at the inner edge at decoupling 
nor the subsequent gas-binary interaction inside this radius are likely to have an appreciable effect 
on the disk evolution on timescales of interest. At the level of accuracy of our idealized models, we 
do not expect either of these factors to affect too strongly the observable properties of the viscously 
spreading disk, to which problem we now turn. 



3. Observable Features of the Time-Dependent Afterglow 

3.1. Bolometric Light Curve 

The power per unit area viscously dissipated from each face of a thin Keplerian disk is equal 
to F ~ (9/8)i^EQ^. Substituting our expressions for u oc r" and il oc r"^/^, and ignoring advective 



speeds are, however, generally subsonic for realistic E profiles, as long as dlnT,/dlnR ^ lO"' at decoupling. Note that 
the step-function profile shown in Figure[na) rapidly evolves to a profile with subsonic radial speeds throughout. 



-13- 



loss, the power radiated by the disk is 

^ 0.059 LEdd X 5>-2A-i X"°-"'(/?-iCr''-''^'' °' f ^^r"-^dr, (3.1) 

J So 

where Lgdd = A-ncGmpM / is the Eddington luminosity for an object with a mass equal to that 
of the binary. The dissipated power depends weakly on a and 9, scales roughly linearly with the 
binary mass, and is most sensitive to the gas distribution parameters S and A. At decoupling, the 
second integral in equation (|3.1I) evaluates to a value ~ 1/(1 -m), and thus the initial luminos- 
ity of the disk at decoupling is expected to be ~ 0.04LHdd for fiducial parameters and m ^ -1/2. 
In integrating equation (13.11 ) before the binary merger, we do not calculate emission from radii 
inside R < 2Xa{t). In this region, the potential is highly time dependent and non-axi symmetric, 
and the assumption that the gas is in nearly circular Keple rian orbits about the binary's cen 



ter of mass bre aks down. Hydrodynamic al simulations (e.g.. iMacFadyen & Milosavljevidl2008 : 



Hayasaki et al.l 120081: iCuadra et al.l 120091) suggest that the gas at these radii will form relatively 
dim, quasi-periodic accretion flows around the individual MBHs. We do not address these circum- 
secondary and circumprimary accretion flows in this pape r . The potential role of circumprimary 



gas inside of the binary's orbit is discussed by IChang et al.l (|2009|) . 



In Figure [3] we present several bolometric rest-frame light curves for the evolving disk, for 
three different values of the viscosity power-law index n, using equation (12.71 ) for the surface 
density profile at decoupling. We show light curves for two different surface density profiles im- 
mediately outside the disk edge at decoupling: a decretion-like disk (m = -1/2: thick lines), and 
an accretion-like disk (m = 0: thin lines). The vertical scale is the brightening Lvisc(0» in units of 
the disk luminosity at decoupling Lq, while the horizontal scale is the time relative to the merger. 
In each case, the dissipated luminosity increases steadily by over an order of magnitud e withi n 
~ 20 yr of the merger. There is significant evolution before and after our re-evaluation of ImPOSI 's 
estimate for the onset for the afterglow, ?em = (3/4)fGw,() ~ 9Mg-^ yr after the merger. 

The sudden brightening that accompanies the merger in Figure [3]is due to the fact that we only 
integrate the circumbinary disk emission outward of the radius R = 2Xa(t), as described above: 
since the binary will merge within ~ 48 hours once it reaches a separation of ~ 50GM/c^, the 
integrated area rises sharply just before the merger. This is a somewhat artificial ingredient of our 
model. We note, however, that at merger the depth of the potential well near the center of mass does 
increase drastically, as does the area where the gas is free to spread viscously. This sudden change 
in the central potential well could be associated with significant emission, for i nstance if even trac e 
amounts of gas remain present between the components of the merging binary (|Chang et al. I I2OO9I) . 



The disk emission can exceed LEdd(^), the Eddington luminosity limit associated with the 
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binary mass M, without violating local Eddington flux limits. This result may appear contradic- 
tory, but it can be understood as follows. The Eddington flux limit associated with the vertical 
gravitational field of the disk (ignoring self-gravity) is 



1 + 



' c GM 



i?2 



(3.2) 



where H is the scale height of the disk and ~ 0.40 cm^g"^ is the electron- scattering opacity. 
The maximum luminosity of a steady disk is then L, ~ 47r J FEdd,z R dR, where we have introduced 
an additional factor of two to account for the fact that the disk emits from two faces. The leading 
geometric factor in equation (13.21 ) is a monotonically increasing function of H /R. It approaches 
the value H/R in the limit H /R <^ 1, is of order ~ 0.7 when H /R ~ 1, and has an asymptotic 
maximum of unity in the limit — )• oo. The luminosity L emitted by an annulus with inner and 
outer radii i?in and i?out is capped by the inequality 

L< r°\rr^R-' dR = \n(^)LEd,. (3.3) 

Thus, the bolometric luminosity can reach several Eddington luminosities for reasonable disk sizes. 

We find that the evolution of the light curve is largely determined by the viscosity profile in 
the disk and has a relatively weak dependence on the surface density profile beyond the disk edge. 
Figure |3] suggests that the accretion luminosity of the viscously spreading inner disk may be a 
viable candidate for observational follow-up almost immediately after the main GW event. If the 
electromagnetic accretion signature associated with the spreading inner disk can be detected in the 
"rise" phase, the light curve may provide clues about the nature of viscosity in the accretion disk. 
We find that our disk solutions can develop bolometric luminosities in excess of the Eddington 
luminosity of the merged MBH remnant. These luminosities can be generated even if the locally 
dissipated flux is sub-Eddington (§ |3.4| below). and are enhanced strongly if the value of the mass- 
accumulation parameter S greater than unity. 



3.2. Spectral Evolution 

We now turn to the spectral evolution of the viscously spreading disk. Below, we present 
an abbreviated derivation of the therm al struc ture of the inner circumbinary disk and discuss its 



spectral features. We refer the reader to lMPOSi and our Appendix |A] for more detailed derivations. 



The spectrum differs from that of an exact blackbody because the photons of different fre- 
quencies are thermalized at different depths above the disk. This thermalizing region is referred to 
the "thermalization photosphere" or the "effective photosphere." The temperature at the bottom of 
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the photosphere is equal to the effective temperature as seen by an observer far above the photo- 
sphere, viz. the effective optical depth above this height is of order unity. We call this temperature 
the photospheric temperature, Tp. 

There are two sources of opacity in the photosphere: electron scattering, which has the 
frequency- indepen dent opacity k^s', and absorption, which is dominated by the bound- free process. 



We follow IMP05I and prescribe a Kramer's functional form for the absorption opacity Kabs,u 
pT-y^MO, where /^(O = ^^(1 -e'^) and ^ = hu/kT. We scale f^ahfi.u SO that its Rosseland mean 
recovers the standard Kramer's bound-free opacity for solar metallicity, 1 .6 x I0^'^(p cm^g"' )(T /K)"'/^ cm^g" 
For a wide range of system parameters, we find that scattering is the dominant source of opacity 
in the photosphere at the inner edge at decoupling. In this limit, n^b^.u ^ f^es, we find Kahs,u ~ 
Kabs,*/i^^^5 with Kabs,* ~ 4.7 X lO^*' cm^g"^ X (Q s)(rp/K)"^^''"^ a frequency-independent quantity. 
(See Appendix |A] for the intermediate calculations for the absorption opacity, including the gen- 
eral case allowing for K^bs.iy ~ ^^es.) 

The so-called "graybody" flux ernitted by each face of the disk is given by the expression 



('e.g.. lRvbicki & Lightmanll 19861 : lBlaesll2004l) 



~ TT^— 5,, (3.4) 

where By{Ty) is the Planck function and t^, = /«abs.!//(«;abs,i/ + ^es) < 1 is the ratio of the absorption 
to the total opacity. Both quantities are evaluated at Tp, the temperature at the bottom of the photo- 
sphere. Since Kabs.i/ is a monotonically decreasing function of u, the photosphere emits increasingly 
less efficiently at higher frequencies, relative to a blackbody with the same temperature. In other 
words, the photosphere has a higher temperature relative to a blackbody with the same radiant flux, 
and thus emits at higher frequencies. The expression given in (13.41) recovers the blackbody flux 
iv = T^Biy in the limit where absorption is the dominant source of opacity, i.e. when ^ I. This 
limit is relevant for our disk solutions, in which can span the full range between and 1 in radii 



and frequencies of interest. Note that the flux expression used by IMP05| is lower than ours by a 
factor of two. 

If again advective losses are ignored, integrating the flux of equation (13.41) over frequency 
must give the flux viscously dissipated by half of the disk, (9/S)i'T,Q^. We obtain 

F ■ 9 

"^'^ / F,di^ = E(Tp,n)aT^' = -uJ:n\ (3.5) 



where 

15 r°° 2eJ^ e<d^ 4 
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is the deviation of the bolometric flux from blackbody and e^(Tp, Q) = Kabs.* / i^ahs,* + i^ea) < 1 ■ The 
above approximation for S is accurate to within 10% in the range < ^ 0.9. (See Appendix lAl 
for a more accurate fit to the integral.) 

To calculate the spectrum, we solve equation (13. 5|) numerically for Tp, substitute the result 
into equation (13.41) and integrate: 

' B,{Tp)rdr. (3.7) 

The temperatures at the midplane and the bottom of the photosphere are related through 

T' = ^r/, (3.8) 

if the region between the two heights is optically thick and can be treated as a one-zone gray 
atmosphere. Here r = Oks^II is the scattering-domi nated o ptical thickness between the two heights, 



an d 6 < 1 is a po rosity correction factor. We adopt ImPOSi 's interpretation of the simulation results 



of lTurned (|2004 ) and use 9 = 0.2 as our fiducial value, and define ^0.2 = 0/0.2. 



Neglecting the advected flux, and in the limit where the circumbinary disk is dominated by 
radiation pressure and electron scattering, we may estimate the frequency for the peak monochro- 
matic flux at decoupling, and at late times after the merger. At decoupling, most of the emis- 
sion comes from the inner edge, and the spectrum in equation (13.71) peaks at the Wien frequency 
^ = hu/kTpfl ^ 2.8, from which we estimate 

/^z/peak ^ 15 eV X a«f 5"-^^A^ iM,-°-^2(/?_iC)-''-'^C- (3-9) 

Long after the merger, the inner disk approaches a quasi-steady accretion track satisfying SnuH = 
A/out(l - a/-/?*/-/?), where Mout is a radially constant mass supply rate of the inner disk. The bolomet- 
ric emission is brightest where the quantity 27ri?z^Sfi^ is maximal in the disk, which corresponds 
to i? = 25/16i?*. Applying this to equation (13.51) . we calculate the photospheric temperature and 
find the corresponding peak frequency for the monochromatic flux at very late times to be 

/ A/f n 1 \ / D 2 \ -18/17 

Ke^-0.71keVxM-/-(-^^) (^) . (3.10) 

Above, the dependence on S is replaced by a dependence on the mass supply rate, as at late times 
the disk loses memory of earlier accumulation near the inner edge. For our choice of = 3GM/c\ 



the radiative accretion efficiency is = 1 - a/1 - 2GM/{R^c^) ^ 0.12 (INovikov & Thomdll973|) . 
The circumbinary disk is expected to reach this quasi-steady track after at least ?visc,o ~ 120Mg^ yr 
has elapsed since merger. 
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FigurelH shows the spectrum of our fiducial model disk around an equal-mass, IO^Mq binary 
as calculated at different times with respect to the merger, from decoupling time (~ 3 yr before 
the merger) up to 120 yr after the merger. We show the exact results using equation (13.71 ) in solid 
line^ , and the blackbody spectrum in dotted lines. In the figure, we only account for emission out 
to a radius of lOOOGM/c^, as we do not expect significant contribution of high-frequency photons 
from regions far outside Rq. We find that the spectrum of the circumbinary disk does not evolve 
significantly from that at decoupling until several hours before the merger, because the deepest 
parts of the central potential are not accessible to the c ircumb inary gas until that time. Our initial 
and final spectra agree qualitatively with those given in lMPOSl As with the bolometric light curve, 
however, significant evolution is ongoing well before and after the nominal estimate for the onset 
of the afterglow. Much of the evolution, however, is complete ~ lOMg yr after the merger. 



3.3. Possible Reprocessing of the X-ray Signature 

The frequencies calculated above suggest that an evolving afterglow at low redshift or low bi- 
nary r nass ma y be observable by existing and future X-ray and UV telescopes shortly after the GW 
event. ImPOSI discussed the possibility that the afterglow could be reprocessed to IR frequencies if 
the merged binary is enshrouded by gas and dust. Here we note that the circumbinary disk itself is 
a plausible candidate to reprocess the UV/X-ray afterglow to lower frequencies. 

Ignoring general relativistic effects , the self -irr adiating flux of the disk with viscously dissi- 



pated flux Fvisc can be written (see, e.g., lBlaesll2004|) 



FirriR) = ^ J 27ri?'sin^ivi,,(i?') dR', (3.11) 

where A is the albedo of the disk and d is the angle between the disk surface and the irradiating 
light ray. In our disk solutions, most of the dissipated disk luminosity Z^isc is produced at small 
radii due to the fact that the viscously dissipated flux is proportional to Q^. Thus, if R is much 
greater than the size of the central emitting region then it is a reasonable approximation to replace 
the integral in the above equation with Lyisc sin??. If one further assumes that A is small and that 
the irradiating flux is large compared to the local value of ivisc, the luminosity reprocessed by 
the outer disk will be Li„ ^ Lyisc sin?9. The light-travel time is short (~ 10~^M^ yr across distances 
~ lO^GM I c^) compared to the timescale for viscous evolution of the inner disk, and thus we expect 
the reprocessing processes below to be effectively instantaneous. 



^ The function peaks at ^ = 3.9, so replace the leading factor 15 eV in equation ( 13.91 ) above with 21 eV to 
derive the locations of the peaks in the plots of FigureH) 
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One way that the factor sin?9 can be large is if there is a warp in the disk that offsets the planes 
of the outer and the innermost regions. The circumbinary disk is expected to be aligned with the 
binary orbit in the binary's vicinity, where tidal torques dominate viscous tor ques. Farther out i n 



the disk, the disk may have a quasi-steady warp with respect to the binary (|Ivanov et al.lll999|) . 
Because the warp would dissipate on timescales longer than the afterglow itself, it could act to 
promptly reprocess the afterglow into longer wavelengths. 

We also find that the evolution of the innermost accretion flow itself could effect a disk ge- 
ometry that is conducive t o reprocess ing. If there is no warp and the emitting region is in the 



equatorial plane, then (e.g. lBlaesll2004|) 



As we show below (Figures [5] and |7]), just outside the regions of brightest emission {R < Rq) and 
around the time of merger our disk solutions are often moderately geometrically thick (H /R>0.1) 
and steeply flared geometrically (large positive dlnH /dlnR). That is, the circumbinary disk itself 
may serve as a shroud of gas that can capture and reprocess the high-frequency accretion signature. 

Either of the two effects considered above could cause the luminous inner disk around the 
binary remnant to irradiate the outer disk at a sufficiently large angle i9 to effect a reprocessed 
IR/optical signature of substantial luminosity. 



3.4. Possible Effects of Advection and Super-Eddington Winds 



Above, we largely confirmed the findings of IMP05L with the additional suggestion that ac- 
cretion afterglows of MBH mergers may be observable somewhat earlier than estimated in that 
study. However, we also find that at late times the spreading accretion flow becomes geometrically 
thick near the MBH remnant, with the scale-height-to-radius ratio H/R formally exceeding unity. 
Here, H = T,/ p = cJQ, where Cs is the sound speed \/^P j p and 7 ^ 4/3 is the adiabatic index 
for a radiation-pressure-dominated gas. The radial profile of p is calculated from the disk temper- 
ature profile, which is obtained through equations (13.51) and (13.81) and the total pressure in the disk 

For some binar y and d isk parameters, the inner edge of the disk may already be geometrically 
thick at decoupling see Table [T]). In such a scenario, the circumbinary gas may be free to 



accrete off-binary plane, and the inner cavity may not be as evacuated as we have assumed here, 
and the difference between the pre- and post-merger disk spectra may be less easily discernible. In 
addition, after the merger, it is also possible for the accreting gas to become radiatively inefficient. 
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as in the advection-dominated accretion flow (ADAF; iNarayan & Yil 1 19941 : lAbramowicz et al 
1988 ) model, which would make it more difficult to observi 



As argued in ImPOSL if horizontal advection is competitive with viscous dissipation then the 
physical disk would b e thinner than suggested b y our simple estimate of the ratio H /R, as in the 
"slim disk" model of lAbramowicz et alj (Il988h . In this paper, we do not attempt to develop a 
self-consistent time-dependent model that incorporates vertical disk structure and radial advection. 
Rather, we proceed below as if the disk evolution can be approximated by the thin-disk evolution 
described by equation (12.61 ) while keeping the above caveats in mind. 



The absolute value of the advected flux can be estimated as (e.g. . lAbramowicz et al.lll988|) 

M P^d 



Qadv - 



p 



dlnp ^dlnT 



dXnR dXnR 



(3.13) 



Substituting equation (|2.4I) for the local accretion rate and equation (IA18I) for the radiated flux, we 
crudely estimate the advected to radiated flux ratio as: 



Q 



adv 



4 P,,,H 



Q 



rad 



3 Si?2^]2 



1 + 



Jln(i/E) 
dlnR 



dlnp dlnT 



dlnR dlnR 



(3.14) 



The simple calculations of the evolving disk spectrum performed in § 13.21 are not hig hly de- 
pendent on the disk temperature (our Figure |4] agrees very well with the spectra shown in ImPOSL 
despite the fact that we have calculated a somewhat lower value of Tq). However, the structure of 
the disk at the inner edge at decoupling and the subsequent evolution of the accretion flow is likely 
to depend much more sensitively on the temperature profile. In the radiation-dominated limit, H 
is proportional to PradS'^Jl"^ so from equation (13.131 ) we see that the advected flux is roughly pro- 
portional to uT^Ti~^Q,~^R~^. This implies that near the inner edge, Qadv c<. T^^^^ approximately (see 
Appendix lAl for how disk quantities such as Q scale with the temperature near the edge). Thus, we 
expect the possible transition of the disk from a thin disk to an advective/slim disk or a geomet- 
rically thick ADAF - and the corresponding spectral ev olution - to be a delicate function of the 
thermal structure of the disk at decoupling. For example. ImPOSI estimated a midplane temperature 
at the inner edge that was ~ 1.3 times higher than we did, with most of the discrepancy due to the 
aforementioned factor of two in the flux formalization. Although this does not result in significant 
differences between our values and theirs for the decoupling radius or the relevant timescales, it 
does lead those authors to estimate that the advective flux is more significant at the disk edge at 
decoupling than we do. We estimate that immediately outside the disk edge, 

Qadv 



Q 



rad 



>lO-^a'Js''X-'''Mf-^\/3.,Cr-'0, 



1.9 ^2.4 
0.2 5 



(3.15) 



^Here, we mean a single-temperature, collisional, ra diatively inefficient flow (as opposed to the coUisionless, op ti- 
caUy thin, two-temperature models in the literature, e.g., Shapiro et aljl976[ Rees et aDll982 : Naravan & Yilll995 ). 
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where iMPOSi estimated Qadv/ Qmd ~ 0.44. The flux ratio is higher at the edge due to the steep 
temperature and density gradients there. The equation above suggests that it is highly sensitive to 
all of the disk and binary parameters save for the binary mass. 

In Figure \5\ we show the radial profiles of several additional disk quantities for the same 
parameters, profile and times as shown in Figure [TJ Along with the surface density, we also plot: 
the simple estimate from equation (13.141 ) for the advected-to-radiated flux ratio, Qadv/ Qmd\ the 
estimated scale-height-to-radius ratio H/R; and F/Fndd, the ratio of the locally dissipated flux to 
^Edd ~ H^^c/ Hq^, the local Eddington flux limit associated with the vertical component of gravity. 
For the last ratio, we find 

Z_ ^ o.06^^Vl + (///i?)-V"->a°f^S'-2A-' ™M6-°°^/3_iC)"°''C'- (3-16) 

Figure[5]confirms that for fiducial parameters, the circumbinary accretion flow remains locally sub- 
Eddington. For fiducial parameter values, we find that a relatively high surface density E/Eq > 
1 (or large S) is required to reach a super- Eddington flux near R = R^, and the disk is thus not 
necessarily likely to develop super-Eddington winds at levels which would significantly affect the 
light curves or spectra presented in Figures [3] and IH At late times, however, the nominal advective 
flux exceeds the locally dissipated flux as the disk becomes geometrically thick, which may lead 
the disk to become radiatively inefficient. 

This hypothetical transition of the innermost accretion flow from a thin disk to a geometri- 
cally thick flow could act to suppress the X-ray afterglow. We estimate very simply the effect 
of advection on the emitted spectrum as follows. We assume that the emitted flux is suppressed 
by a factor / = max{0, 1 - Qadv/ Qiad} < 1, and that the advected energy is not re-emitted. We 
also assume that advection acts to suppress emission evenly across all photon frequencies. The 
resulting advection-limited spectrum is shown in Figure[6] alongside a spectrum calculated without 
accounting for advection (i.e. the same disk spectra as those shown in Figure H]). from decou- 
pling to 20 years after the merger. We compare the two classes of spectra at decoupling, and 1 
month, 2 yr, 5 yr, 9.2 yr and 20 yr after the merger, with thicker lines in the figure denoting later 
times. If the disk were to remain radiatively efficient at all radii and times, the emission at high 
frequencies would steadily increase until the quasi-steady thin-disk track is reached. If instead the 
innermost disk becomes advective, this suppresses the high-energy emission at late times. Our 
results indicate that our fiducial disk model may become increasingly less bright as the inner disk 
becomes increasingly geometrically thick. However, the calculations shown in Figure |6] suggest 
that a brightening afterglow may still be visible in the soft X-ray for a short time ~ 10"^?visc.o after 
the merger. This suggests that the circumbinary gas may have two observable signatures of inter- 
est: an initial brightening phase with a soft X-ray afterglow, and a subsequent dimming phase as it 
transitions into an ADAF. 
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To further investigate the possible consequences of advection on the circumbinary disk evo- 
lution, we turn to a different disk model whose inner edge is already geometrically thick at decou- 
pling. We take 5 = 5 and ^ = 1/3 (binary mass ratio q ^ 1/10), keeping all other parameters the 
same as in the fiducial model. The new choices for S and C, are physically plausible ones. If the 
mass supply rate of the outer disk is comparable to the Eddington limit of the central binary, we 
can expect 5 > 1 at decoupling as a result of mass accumulation. We may also expect ^ ~ 1/10 
mergers to be more commoijzl than near-equal-mass mergers because merging galaxies (dark mat- 
ter halos) have unequal masses and because the mass contrast between the MBHs is expected to be 
higher than that of the host halos (i.e. the observationall y inferred MBH- to-host-halo mass ratio 
relation is a steeper-than-linear function of the halo mass: lFerraresell2002[) . Cosmological merger 
tree calculations for MBH mergers consistently show that contributions to the LISA data strearn 
from MBH coalescences will be dominated b y moderately low-^7 events (e.g., Volonteri et al.ll20()3 : 
Sesana et allboovl : iTanaka & Haimanll2009l) . 



With this new set of parameters, the disk decouples at a time ?gw,o/4 ~ 0.30 yr before the 
merger, when the binary has reached a semimajor axis ao ~ 5AGM/c^. The nominal estimate for 
the onset of the afterglow emission is ?em ~ 0.91 yr after the merger. The inner edge of the disk has 
a somewhat higher temperature. To ^ 8.7 x lO^K. We again use the initial density profile given by 
equation (12.71) . In Figure |7] we show the same disk quantities for this thick disk as shown for the 
fiducial disk in Figure [51 This time, the disk remains geometrically thick and radiatively inefficient 
inside the initial edge radius throughout its evolution. As we did for the fiducial disk model, we 
estimate the advective suppression of the spectrum for the thick disk and show results in Figure [8l 
This simple calculation suggests that a disk that is sufficiently geometrically thick at decoupling 
may not, in principle, exhibit any observable spectral evolution. 

On the other hand, we also find that the viscously dissipated flux of this thick circumbinary 
disk becomes super-Eddington at late times (Figure [Hi). iBegelmanI (|2002[) proposed that accretion 
disks with a super-Eddington flux may be able to stay geometrically thin and radiatively efficient. 
Such accretion flows may also generate powerful outflows and flares. While our simple models 
do not self-consistently treat advection or the super-Eddington regime, our results underscore the 
need for further study of this class of accreting systems. 



Provided that they can overcome the "last-parsec problem.' 
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4. Conclusion 



We have presented a simple semianalytic model for the viscous evolution of a thin circumbi- 
nary disk around a MBH binary, in the final stages of the binary's evolution. Using this model, 
we have estimated the time dependence of the thermal spectrum immediately after the binary 
merger. In what may be the most optimistic scenario, a rapidly evolving soft X-ray signature may 
be observed soon after the GW event, perhaps years earlier than previously estimated. It is worth 
emphasizing that the bolometric luminosity is several orders of magnitude brighter than many of 
the EM counterpart candidates proposed in the literature, and in extreme cases may exceed the 
Eddington luminosity of the central MBH. An important feature of the afterglow mechanism dis- 
cussed in this paper is that it is not strongly limited by the mass of the circumbinary disk. This 



are generally limited by the mass of the gas bound to the recoiling remnant (e.e. jBode & Phinnev 


2007 
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2OO9I; 


Corrales et al. 


2009; 
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Anderson et al. 


2009 
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The reason why an afterglow may be observable so early is because the power per unit area 
dissipated by viscous dissipation scales as F oc uTiQ^. The potential is sufficiently deep close 
to the MBH remnant that in this region even low surface densities can generate more power than 
elsewhere in the disk. A central conclusion of this paper is that enough gas may be able to viscously 
follow the binary to merger. The presence of gas in the vicinity of the binary at merger may have 
significant implications for other proposed EM counterpart mechanisms, as well. Some of this 
gas can reasonably be expected to contribute to circumprimary or circumsecondary accretion disks 
which could generate additional observab l e signatures even if the total gas mas s captured in this 
way is small dArmitage & Nataraianll2005l ; lLodato et al.ll2009l ; Ichang et al.ll2009l) . 



In addition to the possibility that the accretion afterglow proposed by ImPOSI may be observ- 
able earlier than previously estimated, our results raise interesting possibilities about the observable 
features of the post-merger disk at later times. Most of the luminosity is generated in the innermost 
disk (R ^ 20GM/c^), and may be partially reprocessed to IR/optical frequencies by the outer disk 
{R > lOOGM/c^). We have considered two mechanisms that could reasonably reprocess a signif- 
icant fraction of the inner-disk emission: a geometrical warp in the disk just outside the emitting 
region; and the steep geometrical flaring of the vertical thickness of the inner disk as it evolves 
We find that the bolometric luminosity of the time-dependent afterglow can exceed the Eddington 
luminosity of the binary, without violating the local Eddington flux limit. We conclude that the 
EM signature — the unobstructed signature in soft X-rays, as well as the reprocessed signature in 
IR/optical — could become comparable to the Eddington luminosity of the binary. 

Our calculations also suggest that as it spreads, the innermost disk may become geometrically 
thick and thus radiatively inefficient, even if it is neither geometrically thick nor radiatively inef- 
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ficient at decoupling. If the disk is initially geometrically thin at decoupling, then it may emit an 
afterglow that is powered by the newly formed deep central potential before becoming geometri- 
cally thick near the remnant at later times. In such a scenario, an accretion disk around a recently 
merged MBH remnant may provide a unique system where the transition from a radiatively effi- 
cient to an inefficient accretion state can be monitored on humanly tractable timescales. 

In the most pessimistic scenario, it seems possible for the disk to become geometrically thick 
even before the binary merges. In such a situation a disk could behave like an ADAF and it is 
unclear whether a circumbinary cavity should exist at all. Even if a cavity were kept open until 
decoupling, adv ective lo sses may suppress any observable spectral evolution for the disk emission. 



As suggested by lMPOSL horizontal advection coul d also act to make t h e disk thinner than we have 



considered here, as in the "slim disk" models of lAbramowicz et al.l (|1988[ ). A slim disk would 
remain somewhat radiatively efficient and thus could still exhibit an observable evolution of the 
spectrum or luminosity. Another p ossible mechanism to keep the circumbinary disk radiatively 



efficient is super-Eddington fluxes (IBegelmanl l2002h . We find that in the same regions of the 



parameter space where the disk becomes very geometrically thick, the viscously dissipated flux 
also becomes very high. The nominally radiated fluxes near the central regions around the MBH 
remnant can be super-Eddington for physically reasonable parameter values, and this could also 
help the disk produce an observable evolving EM signature. 

Advection-dominated accretio n or a super-Eddington flow may produce powerful outflows 



(e.g.. iBlandford & Begelmanlll999|) . The outright viscous dissipation of super-Eddington fluxes 
near the center could also result in a strong outflow. The spin of the MBH remnant and its orien- 
tation are likely to be well constrained from the preceding LISA observation of the merger GWs. 
As such, the observation of a jet near a recently merged MBH binary would present an unprece- 
dented opportunity to study a MBH-powered jet where the spin of the central engine is precisely 
and independently constrained. 

We conclude that the observation of the EM accretion afterglow of a MBH merger is likely 
to provide a windfall of empirical constraints on the physics of gas accretion onto MBHs. Such 
an observational opportunity would be unprecedented on three points. First, the "initial condition" 
for the structure of the evolving flow may be relatively well characterized thanks to models for 
circumbinary disks and the theoretical understanding of the orbital evolution of the binary leading 
up to the merger. Second, the accretion flow is likely to evolve on humanly tractable timescales. 
Such an accretion signature will act as a probe of the viscosity in the flow, in contrast with a 
steady a-disk in which the emission is independent of the a parameter. Third, the mass, spin and 
orientation of the central MBH will have been independently determined by LISA. The last point 
is important in terms of emission geometry, energetics, and particularly significant if the accretion 
flow fuels an outflow or a jet as it fills the central cavity. 
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In this work, we have used ideaUzed models with a simple radial power-law prescription for 
the gas viscosity. However, in the ct-viscosity model the viscosity is also a function of surface 
density and disk thickness. In light of the large range of E and H/R values found in the regimes of 
interest in our solutions, it would not be surprising if a more realistic viscosity prescription led to 
significantly different results from those presented here. Even in our highly idealized calculations, 
the richness of the physical problem that these systems represent is readily apparent. Given advec- 
tion, radiation-pressure dominance, super-Eddington fluxes, issues of disk stability and geometry, 
and rapid viscous evolution, the circumbinary accretion flows around a merging MBH binary are 
likely to produce intriguing observational signatures that would serve as unprecedented probes of 
fundamental astrophysical processes. The interesting possibilities raised by our simple models un- 
derscore the need for more detailed investigations of this interesting class of objects. 

We thank Zoltan Haiman for insightful conversations and constructive comments on the 
manuscript. We are grateful to the anonymous referee for comments that helped improve this 
paper. This work was supported by NASA ATFP grant NNX08AH35G. 



-25- 



APPENDICES 

A. Properties of the Circumbinary Disk After Decoupling 

We review below the properties of the inner circumbinary disk at the time of decoupling. De- 
coupling is defined as the time when ?gw> the timescale on which the binary orbit shrinks, becomes 
shorter than /3?visc, the timescale on which the central cavity refills due to viscous spreading of 
the disk. The procedure is as follows. First, the decoupling condition tow ~ Z^^visc relates the gas 
properties at the inner edge to the orbital parameters of the binary in terms of the disk midplane 
temperature. Second, the midplane tempe rature is calculated by making simple assumptions about 



the thermal structure of the disk. IMP05 performed these calculations for the inner edge of the 
disk; we follow their approach closely and reproduce their results while pointing out a few differ- 
ences. We also discuss possible ways to set the radial power n of the prescribed viscosity in our 
time-dependent solutions. 



Assumptions underlying the model are detailed in the body of this paper and in IMP05I . No- 
tation for physical quantities is given in our ^ Fiducial parameter values are M = 10^MqM(„ 
( = 4-q/(\ +qf- = 1, a = 0.1a_i, (5 = 0.1/3_i, A = 1, 5 = 1, and 60,2 = 1- Unless otherwise specified, 
the subscript "0" refers to the value at decoupling when applied to binary properties; when applied 
to gas properties, it refers to the value at decoupling time and at the inner edge of the disk. 



Al. Properties at decoupling as functions of midplane temperature 

For circular orbits with semimajor axis a, the timescale of gravitational- wave-driven decay of 
the binary separation is given by: 

_ a _ 5 c^ 4 1 
'^*= J^-T6GW^^^ • ^^^^ 

The viscous time is 

t^i,^ = = ^ ^ 144yrx a" A^/^rrM — , (A2) 

3 ly akT ^ -1 (> y iqogM/c^ J ' ^ ^ 

where T = lO^K is the midplane temperature. 

Applying the decoupling condition ?gw,o ~ /?^visc,o gives the following quantities in terms of 
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Table 1. Inner Circumbinary Disk at Decoupling^ 



Variable 


Factor 


a-i 


S 


A 


Me 


/3-i,C 


^0.2 


flo/(GM/c2) 


126 


-0.34 


-0.24 


0.70 


0.08 


0.42 


-0.08 


^EM (yr) 


9.2 


-1.36 


-0.98 


2.80 


1.32 


1.7,0.7 


-0.34 


So(gcm-2) 


6.2 X 10^ 


-0.68 


0.51 


-0.60 


0.16 


-0.15 


-0.17 


To (10^ K) 


1.3 


0.19 


0.86 


-1.95 


-0.28 


-0.49 


0.30 


Hq/Rq 


0.17 


0.76 


2.43 


-3.80 


-0.12 


-0.95 


1.19 


Pradfl/ PgSLsfi 


430 


1.67 


4.25 


-7.35 


-0.04 


-1.84 


2.17 


Qadv,o/ Qrad,0 


2.0 X 10-2 


1.52 


4.86 


-7.60 


-0.24 


-1.90 


2.38 



^At decoupling, the variable in column 1 equals the factor in column 2 multiplied by the column 
head parameters raised to the powers indicated in columns 3-8. All quantities except for ao are 
evaluated at the inner edge of the circumbinary disk. We derive power-law dependencies identical 
to those found by MP05. Note, however, that we derive a somewhat lower temperature than MP05, 
and estimate that the disk is significantly less geometrically thick than they did. See Appendix |A] 
for detailed calculations. 
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Fig. 1 . — Snapshots of the time-dependent solutions to equation (|2.5I) for the surface density S/Eq 
for fiducial disk parameters. We evolve the solution by applying (12.61 ) with n = 0.4, which is 
roughly consistent with the physical prescription oc T/Q at late times and small radii, when the 
disk approaches a quasi-steady accretion track satisfying T,i/ oc (1 - ^JR^/R). Times in the figure 
are scaled to the merger of an equal-mass, IO^Mq binary. Panel (a) shows the initial surface density 
profiles at decoupling. Panels (b), (c) and (d) show the evolved profiles at 1 year before the merger, 
at merger, and 9.2 years after the merger, respectively. We also evolve a third profile (dashed lines) 
that qualitatively describes a scenario where the binary continues to open a gap even after the 
nominal decoupling condition (see text for details). In panels (a) and (b), the semimajor axis of the 
binary is shown schematically with black circles (the orbit shrinks by ~ 25% between the panels). 
After the merger, we impose a zero-torque condition at i?* = 3GM/c^; surface density profiles 
at merger (Panel c) are truncated at this radius and subsequently evolved with the new boundary 
condition. All three nrofiles have the same niialitative time-denendent behavior, siippestinp that 
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Fig. 2. — Same as Figure [TJ but with n = 11/170 ~ 0.065, which is the value consistent with a 
physical prescription for viscosity, u oc T /Q, just outside the inner edge of the disk. We find that 
our qualitative findings are insensitive to the value of n, as long as the quantity 2-n is not much 
smaller than unity. 




10 20 
time after merger (yr) 



Fig. 3. — The viscously dissipated bolometric disk luminosity for fiducial parameters, in units of 
the Eddington luminosity. The time is relative to merger. The thick lines depict the light curve for 
density profiles that have a decretion power-law m = 91n(z^S) / dlnR = -1/2 outside the inner edge. 
The thin lines are for profiles with m = outside the inner edge. The light curves are somewhat 
sensitive to the value of n = dlnu/dlnR, but have the same qualitative behavior. The light curves 
shown here correspond to disks with viscosity index n = 0.4 (solid lines), n = 1 1/ 170 (dashed lines) 
and n = I (dotted lines). Our model predicts significant evolution before and after ?em ~ 9Mg ^ yr, 
our re-evaluation of the nominal time proposed by ImPOsI for the onset of the viscosity-powered 
afterglow. The sharp increase in luminosity at merger time is a somewhat artificial effect that 
arises from the fact that we account only for the flux dissipated outside the radius 2Xa(t) (see text), 
which shrinks rapidly in the days before merger. Note that the luminosity scales with the mass- 
accumulation parameter 5 as L oc 5'^, and so could be significantly enhanced if the inner disk is 
STpnificantlv more massive than a comnarahle steadv-state Shakiira-Siinvaev disk. 
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Fig. 4. — The thermal disk spectra emitted at (from left to right): decoupling (3.1 years before the 
merger), 1 month after, 1 year after, 2 years after, 5 years after, 9.2 years after, 20 years after, and 
120 years after the merger. The thick, solid lines show the spectra calculated using the "graybody" 
formulation in equation (13.71) . while the thin, dotted lines show the spectra calculated assuming 
the disk instead emits as a perfect blackbody. For any given snapshot, the two sets of lines show 
different spectral distributions of the same bolometric luminosity. For the g ray-bod y case, the 9.2- 
year snapshot corresponds to our re-evaluation of the nominal estimate by ImPOSI for the onset of 
the X-ray afterglow, and 120 years is the viscous time at the disk edge at decoupling. 
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Fig. 5. — For our fiducial disk model: the surface density, scale-height-to-radius ratio H /R; esti- 
mated ratio of the advected to radiated flux Qadv/ Qrad; the ratio of the local flux to the Eddington 
limit iVisc/^Edd; and the flux contribution dL^i^c/dR = luRF^i^c normalized to L/Rq. The disk is ge- 
ometrically thin and advection is largely insignificant until the merger. After the merger, however, 
the inner disk becomes formally geometrically thick and advection-dominated (Panel d). Note that 
in panels (c) and (d) the vertical scale is different from the other two panels. The steep dips in 
Qadv/ Qrad are due to the advected flux changing sign (see equation 0.1411 ). 



-32- 



^ FT 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




Logjoy(Hz) 



Fig. 6. — Emission spectra for the fiducial circumbinary disk. The dashed lines show the same 
spectra as calculated in Figure IH from thinnest to thickest curves, at decoupling and at 1 month, 2 
yr, 5 yr, 9.2 yr and 20 yr after the merger. The solid lines show the spectra at the same times, but 
after subtracting the estimated advected flux. The disk is initially geometrically thin and advection 
is not significant until after the merger. At late times, the inner disk becomes geometrically thick 
and radiatively inefficient, and this could lead to reduced overall X-ray emission. The disk emits 
briefly emits soft X-rays at intermediate times, ~ 2 yr after the merger. 
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(d) 3 yr after merger 
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Fig. 7. — Same as Figure[51 but for a more massive and geometrically thicker disk with 5 = 5 and 
C = 1/3. The disk is geometrically thick and advection-dominated in the inner regions of interest, 
from decoupling to several years after the merger. The local flux also exceeds the Eddington limit 
at late times (Panel d). 
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Fig. 8. — The spectra of a more massive and geometrically thick disk, with 5 = 5 and C = 1/3, at 
decoupling and at 1 week, 3 months, 0.91 yr = ?em and 3 yr after the merger. As in Figure [6l later 
spectra are shown with thicker lines. The dashed lines show spectra ignoring advective effects, 
while the solid lines show spectra under the assumption that advection suppresses the emitted flux. 
In the latter case, the disk is geometrically thick and radiatively inefficient across the entire region 
inside the initial gap radius, and thus no obvious afterglow may be observable (all curves overlap 
except the one 1 week after merger, when our calculations show a brief afterglow that is quickly 
suppressed by advection). 
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the midplane temperature at the inner edge, To e 



136^x^-fa:fA^/'(/?_lC)^/^ (A3) 



?EM = ^/?W,o~13yrxr-6'/'a:;/'A^/X/3!('c'/', (A4) 
I/O ^ 2.0 X 10'' cm^ X ToXa^f A'^/'M6(/3-iC)^/^ , (A5) 
fio ~ 4.5 X 10"^ X Tl[^o['r^\-^^l'^M-^ i^-xO'^'^ ■ (A6) 

We parametrize the surface density via the arbitrary relation So = 5MEdd/(37r//o), where the 
parameter S may be interpreted as the product of the outer mass supply rate in Eddington units 
and the excess in surface density near the inner edge over the value expe cted for a standard 



Shakura-Sunyaev disk, which arises from mass accumulation near the binary dlvanov et alJI 19991 : 



Chang et al.ll2009|) . The latter factor depends on tidal interactions and prior accretion history, and 
in general is expected to exceed unity. We obtain 

So ^ 7.3 X lO^g cm-2 x T'^J^Z^^^ SX-^^^^ (P-iO~^^^ ■ (A7) 

The pressure near the edge is dominated by the radiation pressure Prad = (4-(j/3c)T'^, if T > 
lO^K. We calculate the scale-height-to-radius ratio Hq/Rq and the radiation-to-gas pressure ratio 
^rad,o/^gas,o at the disk edge: 



^0 hPr.±oHo 1 7^1 



rad,0 



P. 



Ro V ^0 ^0^0 Sofio-^o 

0.056 xTo^S-'A^Me^-iC, (A8) 
kToT.0 kToJ^ln'^ 



gas,0 



■Prad,0 



Hm^Ho T/UmpPrad.o ' 

>2 

rad.O 



^gas,0 fcroSgQ^ 

^ 57 X T^'J'al{'s-'X''/'Mlip.,Cr/' . (A9) 

Above, 7 4/3 is the relevant adiabatic index, and p = S/// is the gas density. (Note that in 
the H /R profiles plotted in Figures |5] and |71 we solve for the scale height H in terms of the total 
pressure /'gas + ^'rad and do not assume radiation pressure dominance.) 



A2. Calculating the midplane temperature 



We now turn to calculating the midplane temperature of the inner edge of the disk at de- 
coupling. The thermal spectrum of the inner gas disk differs from that of a blackbody because 
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scattering is a significant source of opacity. In such a medium, photons travel a shorter effec- 
tive path before becoming thermalized and the absorption opacity is effectively enhanced as (e.g., 



Rvbicki & Lightmanll 19861 : lBlaesll2004l) 



K-sfi,u ~ ft^abs.i. W 1 + (A 1 0) 

where n^bsM and ~ 0.40 cm^g"' are the absorption and scattering opacities, respectively. The 
second approximation in equation (lAlOl) is applicable if electron scattering is the dominant source 
of opacity, which is generally the case for our problem in the spectral range of interest (UV fre- 
quencies and higher). 

The disk radiates through a photosphere, with photons of different energies being thermalized 
at different depths. The bottom of the photosphere is defined as the height where the effective 
optical depth Teff^„ ^ PpHpK^ff^iy = 1, where Pp is the density there and Hp is the scale height of the 
photosphere. Characterizing Hp as the ratio of the sound speed at the bottom of the photosphere to 
the orbital angular speed leads to 

Pp-j 4 , (All) 

where Tp is the temperature at the bottom of the photosphere. 

The main source of absorption is the bound-free process, which in general is a complicated 
function of composition and photon frequency due to its dependence on the ionization stat e of the 



gas. For consistency and simplicity, we prescribe the absorption opacity in the same way as lMPOSL 
so that the frequency dependence has the same functional form as the free-free absorption opacity, 
viz. 

^^abs,. ~ «;abs,R PpT^'^^L 1 87/^(0] • (A12) 

Above, ^ = hu/kTp, f„{0 = ^"^(1 -e"^) gives the frequency dependence, and Pp and Tp are in 
cgs. The coefficient «;abs,v is scaled so that its Rosseland mean yields the Kramer's opacity scaling 
at solar metallicity, with /tabs.R ~ 1.6 x 10^*^ cm^g"^ and the factor in square brackets having a 
Rosseland mean of unity. Substituting for the photospheric density using equation (|A11I) gives a 
cubic equation for the absorption-to-scattering opacity ratio K^, = /tabs.i^/'^es at the bottom of the 
photosphere: 

4s. f 1 + ^) ^ uo'-^n%-''^'MO. (A13) 

We find that in most cases, the vast majority of the disk emission is emitted from in frequencies and 
regions that are scattering dominated, i.e. n^b^.u/ /«es ^ 1 • In this limit, (|A13I) may be solved trivially 

1 /2 

to obtain Kabs,i/ ~ /^abs,*/. J where we have defined for convenience the frequency-independent 
quantity ^abs,* ~ 4.7 x 10^° cm^g"' x (Q s)(rp/K)"'^/'^. However, in some of our solutions there is 



-37- 



significant emission from regions with Kabs.i///«es ^ 1- All of the spectra and radial disk profiles in 
our figures were computed by solving (IA13I) generally. 



The angle-integrated emergent spectral flux (e.g.. lBlaesll2004j) is 



F, = 7T^^B,, (A14) 



1 + eu 

where B^, is the Planck function and e^, = Habs,u / {i^abs,u + i^es) = (1 +K~^)~^ < 1 is the frequency- 
dependent ratio of the absorption opacity to the total opacity. Both e,, and 8,^ are evaluated at the 
bottom of the photosphere. Each face of the disk contributes half of the total flux F. Integrating 
equation (IA14I) . we obtain: 



F 
2 



F^du = EaT^, (A15) 



where S is the deviation of the flux from blackbody, 

„ _ 15 f°° lej^ e-^di 



^'Jo l + ey'(0/.(0 



(A16) 



S may be expressed as a function of the frequency-independent quantity e*(rp, Q) = /iabs,*/('tabs,* + 
Kes) < 1. In the limit el^^ <^\,E^ 30n-^d^^ fu'^^^e-^ ^ O.SlSel^^. In the opposite limit 
~ 1 (i.e. very close to blackbody), S ^ 1/[1 -l-(e;;;^ - 1)^^^]. We find that the brightest regions of 
our disk generally span the entire range < < 1 . We find a useful fitting formula for S to be 

- ' (A17) 



l-0.127ey'l + (e;i-l)""' 

The above prescription evaluates S accurately to within ~ 1% for the entire range of e* while 
reproducing the asymptotic behavior at the extreme limits e* ^ 1 and ^ \. We also find that a 

1 /2 

much simpler formula, S ^ (4/5)e* , calculates E accurately to within 10% for < e ^ 0.9. 

The optical depth between the bottom of the photosphere and the midplane is dominated 
by electron scattering, and thus has negligible frequency dependence. Provided that this region 
is optically thick and can be described by a one-zone formalism, the midplane and photosphere 
temperatures are re lated b y = {4/3)T^/t, where r is the optical thickness between the two 



heights. Following 'MP05L we write r = On^Jl, where 6 = 00.2/0-2 < 1 is a porosity correction 
factor. Ignoring horizontal advection, we equate F with the power viscously dissipated per unit 
area in the disk: 

F = 2x S(e.)^ ^ 3207^ ^ 9 ^^^^^ 
3t 15r 4 
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We now solve for To e by substituting equations (IA5I) . (|A6I) and (|A7I) into equation (|A18I) . A 
useful intermediate result is 



V '^es 

which gives Tog oc (zyoSo^oe*^^^To)'"^^^^- That is, the midplane temperature scales with the flux 
scaling to the +14/25th power, i.e. more sensitively than in the standard black body relation of 
T oc fJ/'*. We finally arrive at 

To,, ^ 1.3 X °°s3^3/^°°A-^'^/2X"'^''(/3-iC)"''/'"^(,.f' (A20) 
W e derive an inner-edge temperature that is somewhat lower than the value calculated by 



MP05L primarily because their prescription for the emitted flux is lower than ours by a factor of 
two. (Compare their expression for Fy with our equation (|A14|) . See also our discussion that 
follows equation (13.41).) As noted above, the midplane temperature is rather sensitive to the flux. 



We calculate that IMP05I underestimated the flux through the omission of the factor of two, and 
slightly overestimated the frequency integral. Overall, we calculate a bolometric flux 60% higher 
than they did, and our disk temperature is lower by a factor 1.6'"*/^^ ^1.3. We derive identical 
power-law dependencies on the system parameters as they did. 

The lower temperature suggests a slightly longer timescale for the nominal onset of the X-ray 
afterglow after merger, ^em ~ (3/4)/3fvisc,o T^^^'^ ~ 9 x Mg^ yr for fiducial parameters. We also 
calculate lower values for the quantities Hq/Rq and Pmdfl/PgasXh which are highly sensitive to the 
disk temperature (equations PASII and PAQII ) and thus significantly reduced when the factor of two 
correction to the flux is included. We summarize our results below, and in our Table [T] 

a, ^ 126^xa:l'/'^-4'^/2<'OA'/i«Mf^(/3_iC)"/'°^o."^'°°. (A21) 

Qo ^ 5.1xlO-^->xa!;/i°^i4^/^°«A-^i/2X''^''(/?-iC)"''/'°^of°°, (A22) 
fEM ^ 9.2yrxa:\'/''s-''/''X''/'Mf^''/3T\'^'%T\ (A23) 
So ^ 6.2 X lO^g cm-2 x aZl'^^'S^'^'^'X-'^'M'/'' {p.^Q-"^' 0-',^'', (A24) 

^ 0. 17 X aT's'''/'''X-''/'Mf'' iP-yC)-''^'' el''J''\ (A25) 



Ho 
Ro 

■Prad.O 
-Pgas,0 
.1/2 



430 X a!f^°°5l'"^''/400A-14V20^-V25(^_^^)-147/80 ^867/400^ ^^26) 
0.11 X alf /^''5-23V200;^21/10^1/25(^_^^~,21/40^-2V200_ ^^27) 



We calculate a somewhat shorter scale height for the disk than iMPOSl . who had derived Hq/Rq ~ 
0.46 for fiducial parameters. We show in S3. 41 that the advected flux at decoupling is not likely to 
be significant in the circumbinary disk (for 5=1). 
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A3. Prescribing a value for the viscosity power-law index n 

In our model we prescribe a simple radial power-law u oc i?" for the viscosity. In the cigas-disk 
model, however, v has a physical definition with v oc aPgs^/{pVl) oc T /il. This suggests that a 
reasonable value for n should satisfy the relation 

n^^}^ + 3/2. (A28) 
dlnr 

Below, we apply the above relationship between u and T to evaluate an appropriate value for our 
viscosity power-law index n. 

For fiducial parameters, near the disk edge at decoupling and inside this radius thereafter, 
E ^ (4/5)e!^^ and ~ Kabs,*/i^es oc QTp Applying these approximations to equation (IA18I) 
leads to Tp^^^ oc z/EQ^/^. Since T oc t^^% oc S'/^Jp, we obtain 

49/68 

/8«-18)/17_ (A29) 

One way of characterizing the accretion flow is through the power-law index m = d ln(z/S) /dXwR. 
Writing Jin E /J In 7? = m-n and substituting this into equations (IA28I) and (IA29I ) leads to the solu- 
tion n ^ (30-l-49m)/85. If the disk just outside the inner edge behaves like a steady-state accretion 
solution with m ^ 0, then equation (IA28I) suggests that in this region a reasonable value of n 
h 6/11 0.35. If, however, the gas in this region behaves more like a decretion profile with 
m ^ -1/2, then n = 1 1/ 170 ~ 0.065 would be more appropriate. 

The disk surface density profile outside the edge at decoupling, however, is not particularly 
relevant to the viscous refilling rate of the central cavity. Rather, the pertinent value of n is that 
where the torque gradient is highest, where the most rapid evolution occurs, i.e. near the inner edge 
of the disk. However, in this region d \n{T Ini? > is a rapidly changing function of time and 
radius. Thus, there is no single value for n that can fully describe the viscous evolution of an a-disk. 
At late times {t > ?visc(^) after the merger), however, the inner disk evolves toward a standard quasi- 
steady accretion solution with i/S oc (1 - \JR^/R), where i?* is the inner boundary radius imposed 
by the MBH remnant. In this limit, equation (IA28D leads to n ^ 6/17 -I- 49/(170/ a/^?/^* - 1), 
which evaluates to « ^ 0.4 at radii 7? ^ i?o for 1 < R^c^/GM < 6. We thus select n = 0.4 as our 
fiducial viscosity power-law index, as it is consistent with the asymptotic evolution of the inner 
disk at late times. 

Upon experimenting with several values for n, we have found that the choice of n does not 
qualitatively affect the main conclusions of our study, as long as n ^ 1 . The viscous evolution is 
driven by the single timescale r oc ?visc/(l ~n/2Y, and the bolometric and monochromatic lumi- 
nosities also have weak direct dependencies on n. 



T{r,t)(x 
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B. Green's Function for the Viscous Evolution of the Disk Surface Density 



Equation (12.31) . a second-or der diffusion equation, i s line ar if v d oes not depend on S. Be- 
low, we follow the formalism of Lynden-Bell & Pringld (119741 : see also lOgilviel 120051) to derive a 
Green's function solution for the viscous evolution. 

Suppose the solution has the exponentially decaying form = exp(-AO^''^^(^), where 

<j(i?) is an arbitrary function of radius and p is an arbitrary real number. With the additional 
assumption that v behaves as a power-law, v = uo{R/Rq)", equation (12.31) may be rewritten as a 
modified Bessel equation: 



(p + n){ -^/^~"+P + '' + \ 



^ = 0. 



(Bl) 



Selecting p = n- 1 /4 for convenience and making the substitution A = Suok^R^", the equation above 
has a set of solutions 



-In 



, kR'-"^^\ fkR'-"l^ 

A(k)Ji/(4-2n) Z JT +5(fc)Fi/(4-2„) 



(B2) 



l-n/2j ^ ' \l-n/2^ 

where n <2, J and Y are the ordinary Bessel functions of the first and second kind, respectively, 
> is their mode, and A(k) and B(k) are arbitrary coefficients for each mode. 

We note here that even with a simple prescription for u, the general problem of the viscous 
evolution of a circumbinary disk around a binary is not easily t reated analytically. For a case 
where the Green's function is particularly simple, IPringld (|l99lh showed that one may account 
for the binary torques by applying a zero-torque boundary condition diuT^R^/^)/ dR = at the gap 
opening radius. For binaries of interest in this paper, the gap-ope ning radius it self is a function 
of time, which greatly adds to the "extreme algebraic complexity" (jPringldl 19911) that is generally 
involved with applying such a boundary condition at nonzero radius. Our problem, however, is 
greatly simplified by the decoupling condition, which allows us to approximate the early evolution 
as if the circumbinary disk experiences no torques from the binary. We thus proceed as if the 
potential is due to a single central point mass, and approximate the gas orbits as being circular and 
Keplerian. We discuss the justifications and the modest evolutionary effects associated with this 
model simplification in ^2.2[ 

The boundary condition of interest, then, is that the torques vanish at the center. (At late times, 
it will be necessary to consider the inner boundary condition imposed at finite radius by the central 
MBH remnant. We shall address this point shortly.) Because }'i/(4-2,i)(fcy) diverges as j — > 0, our 
physical solution requires all B{k) to vanish. The general solution is a sum of all possible modes 
over all A; > weighted by a set of coefficients A{k), viz.: 



POO 

T.{R, t) = R-"-'^^ / A(k)Ji /(4-2„)(^) exp [-3i^oRo"kh] dk, 
Jo 



(B3) 
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where we have defined y = R^-"l^ /{\-n/l) for convenience. 



The function A{k) can be evaluated through the use of Hankel transforms. For a given order 
of the Bessel function ^, Hankel transform pairs satisfy: 



= / (j){k)JeikR) k dk 
Jo 

POD 

(f)(k) = / <i>(R)UkR) R dR 
Jo 



(B4) 
(B5) 



For our problem, we may construct the Hankel transform pair 

poo 

$(3;) = R"^'/'E{y,0)= / [k-'A(k)]j,/^^2n)(ky)kdk 

Jo 



from which we obtain 



k-'A(k) = I [7?'«+i/4S(y , 0)] J,/(4-2n)(ky') y' dy\ 
'0 



POO 

A(k) = / R'"^'/'E(y\0)Ji/(^4-2n)(ky')kydy' 
Jo 

S0;',0) J:/^4-2n)(ky')kR"^'^"^^ dR' 



-1 

Jo 



(B6) 



Combining equations (IB3h and (IB6h . we may write the solution E(y, as an integral function 
of the initial condition 0), 



n 



-1 



x7l/(4-2«)(^'Vl/(4-2n)(^3') ^ dR' 

Go(R,R',t)m',t = 0)^, 
with the Green's function 

Go = (l--j R-'"-"R""rJ Ji/^4-2n)iky')Ji/i4-2n)(ky)cxp{-3uoR-o"et) kdk 



(B7) 



j^- 1 /4-nj^/5 /4j^l+n 

6(\-n/2)uot 



'l/(4-2n) 



1 /(4-2«) 



■^l-«/2^/l-«/2^« 

6(l-n/2)2z/o? 

2^1-n/2^/l-n/2 



exp 



exp 



^2-n _|_ ^/2-n 



(B8) 



where Ig is the modified Bessel function of the first kind and order i, and we have substituted the 
dimensionless variables r = R/Rq and r = 8(1 -n/2)^?/?visc.o- 
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In deriving the Green's function above, we have applied the zero-torque boundary condition at 
the origin, instead of at some finite boundary radius R^. At late times and small radii, r ^ r^~", the 
solutions obtained using equation (IB8I) have the radial dependence dln(uTi)/dlnR ^ (n-2)r^""/r, 
viz. they converge to a quasi-steady profile with z/S oc Sr" constant in radius. Thus, E diverges at 
the center at late times if n is positive. 

A physically plausible model should account for the boundary condition imposed by the MBH 
remnant at late times. Below, we approximate the Green's function for the case where the zero- 
torque boundary condition is applied at a finite boundary radius R^, > 0. In this case the boundary 
condition determines the relationship between the A(k) and B{k) coefficients in equation (IB2I) . viz: 



<;u{R)=R-^"Aik) 



J 



1/(4-2m) 



(ky)- 



Jl/(4-2n)(ky*) 
Yl/(4-2n)(ky*) 



1 /(4-2«) 



(ky) 



Our new Green's function is then 



(B9) 



Q _ ^-l/4-n^/5/4+n/2 / j 



'l/(4-2«) 



(ky) yi/(4_2„)(fcy)exp 



X 



^1 /(4-2«) (ky* ) ^1 /{4-2n) (ky) 
Yl /(4-2n) (ky*) J I /{4-2n) (ky) 



k dk. 



Rl-"k^ 



4(l-n/2)^ 



tT 



(BIO) 



The integrand above is greatest where k is small, and decays rapidly as k increases beyond k^ ~ 
4Rq~^t~^ . The boundary condition, however, has a significant effect on the viscous evolution only 
where r ^ r^~" = y^R^'^ and y ^y*. Thus, the boundary effect is appreciable only where k^ <^ 
4y~^ > 4y~^. In the limit of small argument J(:(x)/Y£(x) oc x^"\ and we obtain 



'l/(4-2n) 



(ky,) y. 



1 /(4-2n) 



(ky) 



Yl/(4-2n)(ky*) Jl/(4-2n)(ky) 



l/(2-«) 




(Bll) 



To leading order, then, at late times and small radii the factor in square brackets in equation (IB 101) 
evaluates to 1 - \JR*/R. This factor may be taken outside of the integral, and gives the asymptotic 
behavior z/S oc 1 - ^/R^JR, which is precise ly the standard sol ution for a steady thin accretion disk 
with a no-torque inner boundary at R^ (e.g.. Frank et al.ll2002|) . 



We do not evaluate the exact Green's function solution at early times or at large radii, where 
the effect of the boundary condition is minimal. Instead, we approximate t he function by ensuring 
the co rrect boundary behavior at late times and small radii, as suggested by lLynden-Bell & Pringle 
(11974). Accounting for the fact that the boundary at i?* exists only after the time of merger, which 
evaluates via the decoupling condition to r^erge = 2/3(1 -n/2)^, we approximate our new Green's 



-43- 



function as follows: 



Gir,r',T) = < 1-W — exp 



2-n 



T^{t Tmerge)_ 



Go(r,/,T), 



(B12) 



where Go is the Green's function given in equation (IB8I ) for the case i?^, = 0, = R^^/Rq, and IZ is 
the ramp function, defined as lZ{x < 0) = and lZ{x > 0)=x. 
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